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Over the past years, there have been many efforts towards generating in¬ 
teractions between two optical beams so strong that they could be observed 
at the level of individual photons. Such strong interactions, beyond open¬ 
ing up a new regime in optics, could lead to technologies such as all-optical 
quantum information processing. However, the extreme weakness of photon- 
photon scattering has hindered any attempt to observe such interactions at the 
level of single particles. Here we implement a strong optical nonlinearity using 
electromagnetically-induced transparency and slow light, and directly measure 
the resulting nonlinear phase shift for individual photons. This is done by il¬ 
luminating the sample with a weak classical pulse with as few as 0.5 photons 
per pulse on average, and using post-selection to determine whether a given 
pulse contained (approximately) 0 or 1 photons. We present clear data show¬ 
ing the quantized dependence of a probe beam’s measured phase shift on the 
post-selection result, for a range of input pulse intensities. We believe that this 
represents the first direct measurement of the cross-phase shift due to single 
photons. 

Modern optical physics has revolved principally about two poles: nonlinear optics, where 
rich effects are generated through the interactions of photons with one another, but which 
because of the weakness of those interactions typically manifests itself only for pulses con¬ 
taining billions of photons; and quantum optics, where phenomena such as entanglement 
have been widely studied, but where photon-photon interactions are negligible. In fact, it 
has been well known in quantum optics since the Nobel Prize-winning work of Roy Glauber 
that in the linear regime, the classical and quantum theories of electromagnetism make iden¬ 
tical predictions (except in the presence of nonclassical sources of light - which themselves 
rely on nonlinear effects)It has therefore long been a dream to move into the realm of 
“quantum nonlinear optics,” where sufficiently strong interactions could create a complex 
many-body interacting quantum system in the context of optics. Such interactions could 
enable the production and detection of novel entangled states including few-photon bound 
states^’'^, as well as new architectures for non-demolition measurement of photon number^“^, 
quantum teleportation®, low light level switching®, and quantum logic gates^®. Decades 
ago, important steps in this direction were taken in the realm of cavity QED^^, leading for 
instance to major advances in the study of entanglementand more recent dramatic 
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developments in the context of superconducting qubits^®’^^. Greatly enhanced interactions 
for photons have recently been observed, using electromagnetically induced transparency 
and slow lightmicrostructured fibres^®, atoms in hollow fibres^®, a single atom strongly 
coupled to a micro-resonator cavity^^, atoms or other impurities near tapered fibres^^’^^ or 
bottle resonators^'^ . Most recently, the application of “Rydberg blockades” has enabled 
huge interaction strengths, already leading to the observation of strongly modified quan¬ 
tum statistics^’^^^^°. Sum-frequency generation of two heralded single photons has also been 
demonstrated^^. In parallel, there has been some controversy about the applicability of such 
strong interactions to quantum logic due to fundamental noise limits^^“^'^ and suggestions 
that intermediate-strength nonlinearities might be a powerful approach^^’^®. At the present 
time, the nonlinear phase shift written by a single photon on a probe beam has to our 
knowledge never been reported, and prior experiments have been performed with classical 
pulses with average photon numbers on the order of a few hundred^^ in free space, or as low 
as 16 or 0.1 in a hollow-core fibre hlled with atomic vapour^*^ or a nonlinear photonic-crystal 
fibre^®, respectively. 

Here we demonstrate that by illuminating a sample of atoms with a weak coherent state, 
but post-selecting on subsequent detection of a photon at the far side of the sample, we 
can observe the nonlinear effect of that one additional photon on a probe beam. Similar 
“intensity-held correlations” were previously used to enhance quantum effects in a cavity- 
QED system^®''^^ . Using electromagnetically induced transparency in a gas of laser-cooled 
®^Rb atoms, we report the observation of a nonlinear phase shift of 18 microradians per 
photon, and unequivocably demonstrate the contribution due to individual post-selected 
photons. The sample was illuminated with weak classical pulses with as few as 0.5 photons 
on average, and post-selection allowed us to subsequently determine whether a given pulse 
contained (approximately) 0 or 1 photons. We present clear data showing the quantized 
dependence of the measured phase shift on the post-selection result, for a range of input 
pulse intensities. This represents the hrst direct measurement of the cross-phase shift due 
to a single photon. 

We use an atomic level scheme^® based on electromagnetically-induced transparency^^ 
(EIT), which allows for very strong near-resonant interactions. EIT is a coherent atomic 
effect in which different excitation pathways interfere destructively, eliminating linear ab¬ 
sorption and simultaneously producing sharp dispersive slopes; see the inset of figure 1. In 
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order to establish EIT, two phase-coherent laser helds (probe and coupling) form a A-system, 
addressing a common excited state. When the two-photon resonance condition is satished, 
i.e. their respective detunings are equal, the medium becomes transparent to these helds and 
each experiences a modihed refractive index prohle. The presence of an additional ‘signal’ 
held inside the medium serves to ac-Stark shift one of the ground states of the A-system, 
pulling the probe and coupling helds out of two-photon resonance. As a result, the probe 
held experiences a change in refractive index, acquiring a phase shift that is proportional 
to the slope of this steep dispersive feature. For low signal power, the induced Stark shift 
is linear in signal intensity; the resulting probe phase shift is, therefore, proportional to the 
number of photons in the signal pulse. Measurement of this phase shift provides information 
about the number of signal photons present in the interaction region. 

Figure 1 shows a schematic of the experimental setup as well as the level scheme. A 
cloud of laser-cooled ®®Rb atoms serves as the nonlinear optical medium, with the probe 
and coupling beams passing through orthogonal to each other. Short signal pulses are sent 
counter-propagating to the probe and both beams are collected by highly reflective beam¬ 
splitters after the interaction. A single-photon counting module (SPCM) is used to detect 
the signal pulses while the probe beam is detected by a fast avalanche photodiode (APD), 
operating in the linear regime; see the ‘methods’ section for the details of the probe phase 
measurement. 

We hrst measure the nonlinear cross-phase shift (XPS) for a range of signal pulse energies 
in order to determine the size of the per-photon effect. Figure 2 plots the cross-phase shift 
on the probe held, versus average number of photons per signal pulse. The phase shift grows 
linearly for low photon numbers and a ht yields a slope of 13 ± 1 /rrad per photon. There is 
a saturation at high photon numbers that occurs when the induced ac-Stark shift becomes 
comparable to or larger than the half-width of the EIT window. The lowest energy per pulse 
that we use here corresponds to an average of one photon per pulse, which is the lowest pulse 
energy ever used for cross-phase modulation in free space. 

Although in the measurement described above we have observed the XPS due to signal 
pulses with an average photon number of 1, approximately 40% of the pulses contain no 
photons at all; and about 25% contain multiple photons. To observe the quantized effect 
of individual photons, we send in even weaker pulses - containing only one-half a photon 
on average - and trigger on subsequent detection of a single photon at an SPCM. These 
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so-called ‘click’ events occur only a small fraction of the time, due to the weak signal pulses 
as well as detector inefficiency. Because of this finite efficiency, the absence of a click does 
not preclude the possibility that there had, in fact, been a photon in the interaction region. 
Nevertheless - and neglecting for the moment background counts and multi-photon events - it 
can be shown that the best estimate possible of the photon number in the interaction region 
increases by exactly one when the detector hres, regardless of its efficiency. An intuitive 
picture for this involves recognizing that for an incident coherent state, all the “undetected 
modes” are uncorrelated with the “detected mode”: gaining the information that there is 1 
photon, rather than 0, in the detected mode thus has no effect on how many photons one 
should estimate are present in all other modes. One can therefore observe the effect of one 
additional photon on the probe beam by comparing cases when the SPCM hres with cases 
where it does not. 

When background and multi-photon events are included, one can still calculate the shift 
in inferred average photon number due to a ‘click’ event, and it remains close to (but no 
longer exactly equal to) 1; see the Supplementary Information for details. This difference 
depends on the incident average photon number \a\‘^, the overall collection and detection 
efhciency r], and the background photon detection probability Pb- When the probability of 
clicks due to signal photons is much higher than that of ones due to background photons 
- but still much lower than 1 {Pb 1) - the difference is close to 1. However, at 

low incident photon numbers false clicks due to the background photons cause the inferred 
difference to be smaller than 1. On the other hand, the difference can be larger than 1 when 
the contribution of higher photon numbers to the click events has to be included. Figure 3 
plots the inferred average photon number in the interaction region, riinf, versus the average 
incident photon number, \a\‘^, for click and no-click events. It also shows the photon number 
one would infer from a given number of counts at a photon-number resolving detector, for 
comparison. 

We can separate out the instances in which the SPCM detected a photon (the ‘click’ cases) 
from those in which no click was registered by the detector. This allows us to measure the 
nonlinear phase shift of the probe in the two cases separately, see hgure 4. For 0.5 incident 
signal photons per pulse, the inferred average photon number in the interaction region for no- 
click events is roughly 0.3 whereas it is approximately 1 for click events; see the square shaded 
in green. The phase shift measured for the no-click cases is statistically indistinguishable 
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from zero (2±3 /irad), while the click events result in a non-zero phase shift of —13±6 /irad, 
consistent with the per-photon nonlinear phase shift of —13 ± 1 prad inferred from the £t 
to the data in hgure 2. 

For the remaining data points, the average incident signal photon number and/or the 
center detuning of signal pulses was varied. The most signihcant feature is that the nonlinear 
phase shifts for click events are always larger in absolute size than those for the no-click 
cases. The magnitude of the phase difference between click and no-click events averaged 
over all data points, i.e. the nonlinear phase shift due to a post-selected single photon, is 
— 15 ± 3 prad (—18 ± 4 prad after correcting for hnite background; see the Supplementary 
Information for details). 

In order to conhrm that the observed effect is not due to systematics, we have also taken 
data in the absence of signal pulses or atoms, and with large signal detunings. The data in 
the region shaded in blue displays these checks for systematics. The most important feature 
to highlight is that the click- and no-click probe phase shifts for all systematic checks are 
equal to within the error-bars. 

Finally, the inset of hgure 4 plots all the post-selected data versus inferred average photon 
number corrected for the variable signal detuning; see the Supplementary Information. The 
solid line has a slope of —14 ± 1 /irad per photon which is inferred from the £t in hgure 2, 
the optical density, and the detuning dependence of the nonlinear ehect. We see excellent 
agreement between the data presented here and the value of XPS per photon extracted from 
hgure 2. 

The theoretical prediction for the magnitude of the cross-phase shift, given the param¬ 
eters of our experiment, is 13 prad based on the model presented in'^^; see Supplementary 
Information for more details. The values we measured in our experiment are in very good 
agreement to the theoretical value. 

In summary, we have used FIT to make the first observation of cross-phase shift due to 
free-space signal pulses containing 1 or fewer photons on average. Moreover, by conditioning 
on detection of a photon in the signal pulse after it is transmitted, we have been able to 
detect the (quantized, in the ideal case) ehect of a single additional photon on a probe beam, 
and measured this single-photon phase shift to be —18 ±4 prad, consistent with the results 
for classical signals. This is a step towards the development of further ’’quantum nonlinear 
optics” techniques, including the potential for scalable quantum logic gates. 
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Figure captions 

Figure 1. Schematic of the experimental setup. Counter-propagating probe and signal 
beams are focused to a waist of 13±l/im inside a cloud of laser cooled ®^Rb atoms conhned in 
a magneto-optical trap. High reflectivity beam-splitters (90% reflectivity) are used to collect 
the signal and probe beams after the interaction. A collimated coupling beam, propagating 
perpendicular to both probe and signal beams, creates a 2 MHz EIT window for the probe 
(see the inset). The amplitude and phase of the probe are measured using frequency-domain 
interferometry. The single-photon detection events on the SPCM are registered as tags in 
the amplitude of the probe. The inset shows the level scheme used. 

Figure 2. XPS versus average photon number per pulse. The nonlinear phase shift 
depends linearly on the photon number at lower intensities. A £t to the low-photon-number 
data yields a slope of 13 ± 1 prad per photon while the deviation at higher photon numbers 
arises due to higher-order nonlinearities. The inset shows a typical linear phase prohle 
(green) and optical density (red) as seen by the probe with the arrow indicating where the 
on-resonance component of the probe laser is locked. Other relevant parameters include 
signal center detuning = -10 MHz, OD = 2, EIT widow width = 2 MHz. 

Figure 3. Inferred {riinf) versus average photon number in the interaction region. The 
overall collection efficiency is assumed to be 20% and the background click rate is taken to 
be 10% for solid and dotted lines. The circles show the photon number values inferred for 
the data points in figure 4 for no-click (red) and click (blue) events. The overall efficiency 
percentage for each data point (numbers beside circles) is slightly different which accounts 
for the discrepancies between the data points and the solid curves. The average photon 
number in the interaction region for the data points is lower than the incident photon 
number because of the finite signal absorption. The dotted green lines show the photon 
number which would be inferred were a number-resolving detector used. The solid blue line 
could also be obtained from a weighted average of the dotted lines with non-zero number of 
clicks. 

Figure 4- Post-selected single-photon XPS. Most notably, for an average incident photon 
number of 0.5 (green-shaded region), the XPS for no-click and click events are 2 ± 3 and 
—13 ± 6 prad, respectively, which definitively shows the effect of a single post-selected 
photon. For the other data points, the average incident photon number and/or the signal 
center detuning is varied. Taking all the data points together, the magnitude of the post- 


7 


selected single-photon XPS is —18 ± 4/irad. The inset shows the post-selected XPS versus 
Uinfi^TT X 18 MHz)/As, inferred average photon number corrected for the variable signal 
detuning. The solid line has a slope of —14 ± 1 /irad per photon. Other relevant parameters 
include EIT window = 2 MHz and OD = 3. The data in the region shaded in blue are tests 
for systematics as explained in the text. 


Methods 

Atom preparation. The atoms are prepared in a magneto-optical trap using three beams 
and their retro-reflections along with a magnetic held gradient of 20 G/cm on the axis of 
the quadrupole coil. Each of the three beams has a ‘trapping’ (2.5 cm diameter) and a 
‘repumper’ (1 cm diameter) component. The trapping beam is tuned to 20 MHz below 
the cycling transition F = 3 —?■ F' = 4 and the repumper is tuned close to resonance on 
F = 2 ^ F' = 3 transition to keep the population in the trapping ground state. The 
beams and the magnetic held gradient are turned off every 22 ms, leaving the atoms to 
freely expand and allowing them to be probed for 1.5ms, as explained in figure 7. 

Probe phase measurement. The nonlinear effect we are interested in is a change in the 
probe refractive index, which is linear in the signal photon number and manifests itself as 
a phase shift of the probe held. In order to measure this probe phase shift, frequency- 
domain interferometry is used. The probe beam is comprised of two different frequency 
components, which copropagate through the atomic cloud. One frequency component is 
tuned on resonance with the F = 2 —)■ F' = 3 transition while the other is -1-100 MHz 
detuned, serving as a phase reference; any phase or amplitude change due to the atoms 
appears as a phase or amplitude change of the resulting 100 MHz beating signal. The probe 
beam is detected on a fast avalanche photodiode; demodulation of the resultant electrical 
signal at 100 MHz allows us to extract the phase and amplitude change experienced by the 
on-resonance probe held. The analysis bandwidth for the demodulation is 2 MHz, which is 
matched to the EIT window for optimum signal-to-noise ratio. A stable 10 MHz clock is used 
to generate a 100 MHz signal used as the reference for the modulation and demodulation. In 
order to eliminate the effects of any slow phase drifts, we use the average probe phase over 200 
ns durations before and after the expected XPS as a phase background, which is subtracted 
off from the XPS. An important advantage of using frequency-domain interferometry and 
demodulation is that the measurement is insensitive to any source of variation that does 
not have a component at 100 MHz. We measure a single-shot phase uncertainty of 50 mrad 
(there are approximately 4500 probe photons per measurement window, which is equivalent 
to a shot noise of 15 mrad). 

Probe and coupling fields. Our master laser is locked 30MHz below the F = 2 —)■ F' = 3 
transition of the D2 line in ®^Rb. A portion of this light is frequency-shifted using an acousto¬ 
optic modulator (AOM) driven at -1-130 MHz to produce the off-resonance component of 
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the probe, which is then passed through a second AOM driven at -100 MHz to create the 
on-resonance component. These two components are then combined on a beam splitter, 
forming the probe beat signal. Another portion of the master laser is modulated using an 
electro-optic modulator (EOM) driven at 3 GHz. This modulated light is used as the seed 
for injection-locking a diode laser; locking to the lower side-band of the modulation produces 
a coupling beam that is phase-locked to the probe, a necessary condition for establishing 
EIT. This coupling beam passes through a third AOM, which serves to frequency shift it 
onto resonance with the F = 3 ^ F' = 3 transition, as well as to shutter the coupling light 
on and off in sync with the atomic duty cycle. The width of the EIT feature is linear in the 
coupling intensity and is measured to be 2MHz for all the data taken for this paper. We use 
TT and cr+ polarizations for the coupling and probe fields, respectively. 

Signal pulses. A portion of the injection-locked laser light is split off to produce the signal 
pulses. An AOM is used to frequency shift the light close to the F = 3 —?■ F' = 4 transition 
as well as to create short pulses (40ns or 100ns FWHM). The temporal profiles of the signal 
pulses are measured on an avalanche photodiode. By integrating the pulse power over time, 
the average energy, and therefore the average photon number, of each signal pulse can be 
measured. The polarization of the signal light is a~^. 

Single-photon detection. The signal pulses are collected into a multimode fiber and are 
detected on an SPCM. Other than the signal photons, any stray light which leaks into the 
signal collection can result in clicks. We use gate signals with durations matching the signal 
pulses, to reject background photons. The residual probability of background counts is 13% 
for 100ns pulses and 6% for 40ns signal pulses. Most of these result from the near-resonant 
scattering of the probe and coupling beams from the atoms. The overall collection and 
detection efficiency is around 20%, and for higher incident photon numbers more attenuation 
is added to keep the total click rate around 20% —30%. The details of the tagging procedure 
are explained in the Supplementary Information. 

In order to reduce the background photon rate for cases of incident average photon number 
of 0.5 and 1, the signal pulse duration was chosen to be 40 ns, while it had been 100 ns 
for the earlier data taken for higher signal power. One might expect this change to make 
the signal pulses more intense and therefore make the nonlinear phase shift more than twice 
as big. However, we have shown theoretically"^^ that because of the bandwidth mismatch 
between the EIT window (2 MHz here) and the signal pulses, the peak phase shift nearly 
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saturates and we expect an enhancement of only 1.5. Therefore, given the size of our phase 
measurement uncertainty, the two values are expected not to be statistically different. 

The inferred average photon number difference between click and no-click cases is slightly 
different from unity due to the effects of background counts, multiple-photon events and 
detection efficiency. In order to calculate the nonlinear phase shift per post-selected single 
photon, the measured phase differences are divided by the values of the inferred average 
photon number difference. 

Data collection. For each data point, we took approximately 300 million shots over 14 
hours, 90 million of which resulted in clicks at the SPCM. Because of our tagging procedure 
90 million shots were discarded, and out of the remaining shots we observed 60 million click 
events and 150 million no-click events. 

Focus size. The waists of the probe and signal beams inside the cloud are 13±1 /rm, 
corresponding to a (two-sided) Rayleigh range of roughly 1.4 mm. The choice of this focus 
size is to ensure that the Rayleigh range of the beams matches the size of the cloud. Focusing 
the beams tighter than the size of the cloud would be detrimental because it would reduce 
the interaction length and increase the probe phase shot noise (the intensity of the probe is 
limited and fixed by the saturation intensity of atoms). Also, focusing the beams any less 
tightly would produce a smaller intensity for given signal pulse energy, thereby decreases 
the size of the nonlinear effect. 
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Appendix A: Supplementary information 


Inferred average photon number in the interaction region. In order to calculate the most 
probable number of signal photons interacting with the probe, one can use Bayesian inference 
to obtain the conditional probability of having had Uph photons in the interaction region 
given the occurrence (“yes”) or lack (“no”) of a signal photon detection event (‘click’) after 
the interaction: 


1 I I ^TLpf^ 

P{nph\no) = —P{no\nph)Pinph) = - -(1 - 

No Uph'. 

1 I I 

P(np,,|yes) = —P(yes|np/,)P(np/,) = ^ ^ (1 - (1 - (Al) 

where Nq and Ni are normalization factors, rj is the overall collection and detection efficiency, 
and P{nph) = exp(—|ap)|ap"'p'*/np/i! is the incident signal photon number distribution with 
average \oi\^. Using the conditional probabilities given above, one can calculate the average 
photon number in the interaction region in each case. 


'^inf,no |Q:| (1 ^)) 

^inf,yes I®! P h 

’^m/,no P 




1 _ 

?7|ap 


Psigiyes)' 


(A2) 


where Pgigiyes) = 1 — Pgig^no) = 1 — is the probability of a detection event. For 

low count rates, 1, the expression above reduces to hj„/^yes ~ Hinf,no P 1; fhe 

difference in the inferred average photon number between click and no-click events is unity, 
independent of both the detection efficiency and the average incident photon number. It is 
straightforward to include the effect of background photons: Psig(no) —)■ Psig(jio)Pbkg(jio) 
and Psig(yes) —?■ 1 — Psig{no)Pbkg(jio) where Pbfcg(no) is the probability of getting no clicks 
from background photons. 

Detuning-dependence of the cross-phase shift. Figure 5 shows the level scheme used 
for this experiment; the presence of an off-resonant signal pulse introduces an ac-Stark 
shift, which effectively detunes the probe held from resonance causing it to acquire a phase 
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shift. Here, we are interested in the dependence of this phase shift on the detuning of the 
signal pulse. So long as the probe samples the linear portion of its refractive index prohle, 
the detuning-dependence of the nonlinear phase shift will be the same as the detuning 
dependence of the ac-Stark shift itself. In the limit of weak signal pulses and weak probe 
held, the ac-Stark shift as a function of signal detuning. A*, is given by. 


A ACS /^2 + (r/2)2’ 

where F = 27r x 6 MHz is the excited state linewidth and is proportional to the intensity 
(and therefore the photon number) of the signal pulse. This detuning dependence has a 
dispersion-like shape that goes to zero on resonance and reaches extrema at one-half of the 
line-width, F. 

A second contribution to the detuning-dependence of the nonlinear phase shift arises from 
the non-vanishing atomic population in the ground state addressed by the coupling held. 
The population in each ground state of an EIT system is determined by the ratio of the 
coupling and probe held intensities. In our case, the probe intensity is sufficient to bring 
some population into the ground state addressed by the coupling and signal helds, which 
leads to hnite signal absorption; see hgure 5. This leads to an ehective number of signal 
photons in the interaction region. 


^eff - 


No 


\ _ g fis(As) 

ds{As) 


(A4) 


where ds{As) = doT'^/ (dA^-i-F^) is the optical density of the signal transition at the detuning 
A^ and Nq is the incident photon number. Here do is the on-resonance optical density. The 
nonlinear phase shift is then proportional to the product of these two contributions. 


A,F/2 

'”A2 + (r/2)2 d,(A,) 

where (j)m is a proportionality constant. 

Figure 6 shows the measured nonlinear phase shift as a function of signal detuning, along 
with hts based on all these contributions. The detuning dependence is measured for two 
different probe intensities, keeping the coupling intensity hxed. The £t parameters are 
do = 4 ± 2, (j)m = 500 ± 100 prad for low probe power, and do = 5 ±2, (^rn = 300 ± 40 prad 
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for high probe power. It can be seen that in the case of higher probe power the nonlinear 
phase shift is smaller because of the higher signal absorption. 

Theoretical value of the cross-phase shift. Previously, we showed theoretically that the 
temporal profile of the XPS expected for a single-photon Gaussian signal pulse with an rms 
duration of interacting with an EIT medium with a response time of r is given by^^ 


4>{t) = exp(-t/r) erf(t/\/2rs - Ts/\/2t)^ 

where erf(a;) = 2/y^ dx' exp(—x'^) is the error function, and 


(A6) 


00 — 


d 


(A7) 


-4A^ TTwg Aeit 

is the integrated XPS per signal photon. Here, P is the excited state linewidth, is the 
signal detuning, a at is the atomic cross section as for the signal, Wq is the beam waist, d is 
the change in the optical density for the on-resonance probe without and with a resonant 
coupling beam, and Aeit is the full-width at half-maximum of the transparency. For our 
experimental parameters (P = 27r x 6 MHz, |As| = 27r x 18 MHz, A/aat = 3000, d = 2, 
Aeit = 27r x 2 MHz, = 40 ns, and r = 250 ns) the XPS has a peak equal to 13 prad. 

The maximum achievable cross-phase shift per photon in the N-scheme, due to group 
velocity mismatch issues, is 


“ 4|A,| A ’ 

as shown by Harris and Hau'^^. For the parameters of our experiment this value is 28 /irad. 
The phase shift we measure is lower than this value because our optical density is not high 
enough to saturate the limit posed by group velocity mismatch. 

Measurement and atom cycle. The atoms are first captured for 20 ms in the F = 3 
ground state using a magneto-optical trap (MOT) with a trapping beam red-detuned by 
20MHz from the F = 3 —?■ F' = 4 cycling transition of the D2 line in ®5Rb. Then, for 0.5 
ms, this trapping light is tuned closer to the F = 3 —)■ F' = 3 transition so as to prepare 
the population in the F = 2 ground state. During this time, the magnetic field gradient is 
turned off (and kept off until the next recapture period) in order to avoid spatially varying 
Zeeman shifts, which would lead to dephasing of the FIT system. After this population 
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preparation stage, all MOT beams are turned off, leaving the cloud to freely expand while 
the probe and coupling fields are turned on for 1.5 ms. The amplitude and phase of the 
probe held is continuously measured during this period. Short pulses of signal light (40 ns or 
100 ns FWHM duration) are sent in every 2.4 /is, each pulse constituting one ‘shot’. After 
1.5 ms, the probe, coupling and signal helds are turned off and the recapturing period begins 
again. The cross-phase shift that the signal pulses write on the probe held is obtained by 
splicing the 1.5 ms trace of probe phase data into 2.4 ps shots and averaging over many 
runs; see hgure 7. 

Tagging procedure. In order to see the ehect of a post-selected single photon, we need 
to select out the shots that lead to the detection of a signal photon after the interaction. 
This is achieved by sending a bright hash of light into the probe detector conditioned on the 
click from the single-photon counting module (SPCM) that is used to detect these signal 
photons. This hash of light appears as a spike in the amplitude of the probe held and allows 
us to pick out the shots that lead to a single photon detection event. 

In order to exclude false positives arising from background counts of the SPCM, we use 
time-gating to accept only those clicks which occurred during the signal pulses themselves. 
A logical AND gate is used to do the time-gating, the output of which triggers the hash of 
light that enters the probe detector. 

Prior to the interaction, signal pulses pass through a polarizing beam splitter (PBS). 
One output of the PBS is further attenuated and sent to the atom cloud, serving as the 
single-photon-level signal pulses described above. The other output of the PBS is directed 
to a fast photo-diode, the electrical output of which goes to one input of the AND gate. The 
TTL output from the SPCM serves as the second input to the AND gate and in this way 
we can exclude clicks that arise from background or dark counts of the SPCM. 

Recalling that our probe detection scheme involves demodulation of the probe signal at 
100 MHz, the flash of light used for tagging must contain a 100 MHz frequency component 
in order to make it through the demodulation. This is achieved by using the output of the 
logical AND gate to switch on an AOM, which is driven at 100 MHz. Picking off the first 
and second orders of light passing through this AOM produces a 100 MHz beat signal that 
is then sent into the probe detector. The AOM driver is switched on for 200 ns, producing 
short bursts of light, which show up as spikes in the probe held amplitude. The output 
of the logical AND gate is delayed so that these spikes arrive at the probe detector during 
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the subsequent shot. For this reason, we exclude all shots that follow a successful detection 
event. An important technical advantage of this method to select out shots containing single 
photons is that it is electrically decoupled from the phase measurement electronics to avoid 
any cross-talk between the two systems. 
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Extended data 

Figure 5. Level scheme. The level scheme used to observe cross-phase modulation using 
the D2 line of atoms. The ac-Stark shift due to the signal pulses, pulls the probe out 
of the EIT condition and this appears as a refractive index change proportional to the signal 
intensity. 

Figure 6. XPS versus signal detuning. The nonlinear phase shift is caused by the ac-Stark 
shift due to the signal pulses. Therefore, it has the same dependence on signal detuning 
as the ac-Stark shift. This scaling also depends on probe power because more probe power 
results in a larger population in F = 3 ground state (see the inset of £g. 1) which means a 
larger signal absorption. The overall effect is broadening and smearing of the dispersion-like 
scaling at higher probe powers. 

Figure 7. Measurement and atom cycle. The ®®Rb atoms are prepared in a magneto¬ 
optical trap using a cycle as explained in the text. During the free expansion time, signal 
pulses are sent into the interaction region every 2.4 /is, one ‘shot’, and the phase of the 
probe is monitored continuously. The interesting quantity here is the difference in the phase 
of the probe for cases that there is a signal single-photon detection, ‘click’, compared to the 
no-click cases. 
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FIG. 1. Schematic of the experimental setup. 
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FIG. 2. XPS versus average photon number per pulse. 
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FIG. 3. Inferred {riinf) versus average photon number in the interaction region. 
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FIG. 5. Level scheme. 
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FIG. 7. Measurement and atom cycle. 
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